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The object of the present paper is to investigate the expressions which present 
themselves in the solution of the problem of the conduction of heat in an ellipsoid of 
revolution. For although the question of the stationary temperature of ellipsoids in 
general has been completely solved by means of the functions introduced by Green 
and Lam^ the corresponding problem of conduction has not been so successfully dealt 
with. M. Mathieu, indeed, in his ' Cours de Physique Mathematique/ has shown 
how to reduce the solution to ordinary differential equations, and for the special case 
of an ellipsoid of revolution has shown how to approximate to their solutions. His 
method, which is novel and remarkable, enables him to calculate a few terms of the 
expressions, but does not afford a view of their general constitution and properties. 
In the present paper the question is treated in a more direct and general manner. 

Choosing with M. Mathieu, as coordinates of a point, the azimuth <j> of the 
meridional section through it and the parameters a and ft of the ellipsoid and hyper- 
boloid confocal to the surface which intersect in the point, it is first shown how to 
transform the general equation of conduction to these coordinates. This equation is 
then satisfied by a series of terms of the form e" xnt cos m$ S-J((3)f2j(a), in which k is 
determined by an equation whose roots are infinite in number. 

The function $ is expanded in what Mr. Todhunter, translating Heine's term, 
calls associated functions of cos /3, and we shall also follow Heine in denoting by 

m 

IV*(fO the expression (ft a — l)*^*""*— &c). The language of harmonic analysis has 
been greatly benefited by Professor Maxwell's introduction of the words type and 
degree into the specification of a tesseral harmonic, though we prefer to replace the 
term type by order. We shall therefore call the product P^(/x) cos m(j> the tesseral 
harmonic of the m i]l order and n th degree, and the factor P,/(/x) the associated function 
of the m th order and n th degree. The expansion of the function & will then consist of a 
series of associated functions of the m th order. 

It is shown that the roots of the equation in k fall into two classes, and that the 
corresponding expressions for £ take different forms, for one of which the difference 
between the degree and order of the associated functions involved is an even number 
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and for the other odd. The values of /2(a) divide themselves in like manner into two 
classes. These values are expanded, in the first instance, in terms of the minor axis 
of the confocal ellipsoid, and afterwards in terms of the major axis, the former series 
proceeding by functions which satisfy the differential equation 

dhb 2 du I n.n + l\ A 



dx 9, x dx \ X s 

Of the two solutions of this equation, which are both finite in form, one S„ is finite 
when x=0, while the other T n becomes infinite. 

The expression B n plays for spheres the same part that Bessei/s function plays for 
circular cylinders, and as Heine has employed the term cylinder-function for the 
latter, it would seem consistent with analogy to use the term spherical-function for 
S^ When fl is expressed in terms of the major axis of the confocal ellipsoid two 
expansions are given, one in terms of spherical functions and the other in terms of 
associated functions, and it is shown that each of these series possesses special 
advantages in relation to particular points which arise in the problem of the con- 
duction of heat. 

The properties of these functions are afterwards further considered. 

The above expansions being found, I have next discussed the system of equations 
which determine h y following more or less closely the analysis which Heine has given 
of a similar system, and it is proved that the values of h are all real and definite in 
position, and that for these values the expansions of 3- and fl converge rapidly when 
a sufficiently large number of terms is taken. 

I then show how to express, by successive approximation, the roots of the equation 
in k in powers of a quantity e which depends on the eccentricity, and have entered 
with some fulness of detail into the numerical calculation of a few of the smaller roots 
and of the corresponding coefficients of the functions & and fl, more especially in the 
case of the first of the two classes into which they fall. Besides these particular 
values, however, the general formulae are given, which will furnish them to a certain 
degree of approximation for all values of m and for any value of h 

With regard to the special problem of the conduction of heat, the boundary con- 
dition is supposed to be either that the surface is kept at a constant temperature, or 
that the body is cooling by radiation. The former is mathematically the simpler, and 
we might imagine it realised in the case of a body kept in the midst of an infinite 
fluid after a sufficient time has elapsed for the surface to take the temperature of the 
fluid. With this assumption the different values of X might be found from the equa- 
tion I2 m *=0; and the roots of t T2 °=0are investigated up to e 4 . The general condition 
of radiation is next considered, and it is shown how it may be brought theoretically 
within the range of analysis. I have not, however, thought it necessary to do more 
in this direction than indicate how the successive approximations may be found. 

Although the calculations were undertaken in the first instance for the case of an 
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ellipsoid whose major axis is that of revolution, it is proved that the expressions and 
results may be easily Warred to the ease ef a planetary ellipsoid. And „ the ehief 

obiect of the paper is to examine the forms and properties of the functions which 
present themselves in the mathematical treatment of the subject, I have not thought 
it necessary to enter on the discussion of the movement of heat in a shell bounded by 
two confocal ellipsoids, or in an infinite solid from which an ellipsoidal cavity has been 
taken—questions, however, which might be easily treated by introducing the functions 
T n as well as S # . The analysis would seem also to be applicable to other physical 
preblenr, relating te ellipse ef revolution. 

Perhaps I should explain that if I have entered somewhat minutely into points 
which are not new, I have done so for the purpose of rendering the argument clearer 
and more coherent. 

2. When a solid body is heated in any manner and left to cool, the equations which 
have to be satisfied are, first, the general equation of conduction 

cPY <PV d*V IdY 

* ■— {•*• mJmm — — — """""" — ill 

dx*^ dtf* d?} t dt v 1 ' 

wherein V is the temperature at any point, I the "thermal diffusivity" of the body; 
2° the boundary condition, which we shall suppose to be one or other of the forms 

V=0or^+1jV=0 ......... (2) 

dY dY dY dY 
wherein — =J— +*m— +flnr» I m n being the direction cosines of the normal N to 

dE dx dy *" ° 



the surface measured outwards, and ^ is a constant. 

The first of these conditions corresponds to the case where the surface of the body 
is maintained everywhere at uniform temperature, for we may suppose the zero of 
temperature so chosen as to coincide with the given temperature, whatever it may be. 
The second is the usual condition of radiation according to Newton's law of cooling. 

We have also a third condition, that which gives the initial state of the body : — 

V=V ==/(#, y, z), when £=0 . . . . . . . . (3) 

The general course of solution, whatever the solid may be, is to put 

Y=%A.e~* K \v ....*.♦.. * (4) 
where v is a function of x y z, so chosen as to satisfy the equation 

dx^dy^dz^ A ° ......... {0) 
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and the boundary condition 



dv 
v=0 or — +t)^=0 (6) 



When the appropriate functions v have been found, in general triply infinite in 
number, to satisfy (5) and (6), the constants A may be determined from condition (3). 
Now it is obvious from the nature of the case that any solution of the equations which 
satisfies (1) and (2) and reproduces (3) when £=0, will be the solution sought, and the 
same conclusion can be readily demonstrated by analysis. With regard to equation 
(6), it serves two purposes : first in enabling us to select the appropriate form of v, and 
secondly in furnishing the values of X, which determine the types of heat-movement 
which take place. Poisson has shown, in a very elegant way, that the values of X are 
always real, and as his results are of importance as showing also how the constants A 
are to be found, I shall here reproduce them. 

Let v and v be two functions of x y z satisfying the equations (5) and (6), and let 

d?_ jP d? 
dx 2 dy 2 dz* 



V 2 stand for T^+TTa+y* th en 



(X' 2 — X 2 ) w'dxdydz = Uv V H — v V H')dxdydz 



[' 



, dv dv f , , 
V m~ V dX ]d8 > 



ds being an element of the surface ; hence, if either form of (6) be true, 



\w'dxdydz=0 (7) 

It follows from this equation that the equation in X cannot have imaginary roots; 
for, v being always a function of X, if X=p+?v / "~l^ there will be another root 
\'z=.p— q^/— 1, and the corresponding values of v will be respectively L+M^/— -1 and 
L— Mv/^T. Equation (7) now becomes 

\(U+M*)dxdydz=0, 

which is clearly impossible. We may also employ (7) to find A, for 

Y =$A.v 

A= M~ , dE=dxdydz (8) 

in which the integration is extended throughout the whole of the solid. 



OF HEAT IN ELLIPSOIDS OF BE VOLUTION. 121 

3. For the case of a solid sphere whose radius r=r , these equations become, in 
polar coordinates, 



cPY 2 dV 1 

9.+;: ^ -r 



^ r 2 r ^ 7 



,2 






1 rfV 



V=0 or -T+ ljV=0, when r=r , 

tt/i 
V =/(n /*, <£), 

c£E= —dr.r 2 d(f)diJL, where /x= cos 0, 

The general type of the solution is 

Y=e~ k * it (A cos m</)+B sin m^)P,/B ; , (9) 

in which m may have all integral values from up to co, P,/ cos m</> is the tesseral 
harmonic of the m th order and n th degree, and B^ satisfies the equation 



u » 



+;t+(x s -^)R.=o (10) 



The two particular solutions of this equation are 

,~ /l d\ M sm\r rn /l ^ \ w cos \r 

(as will be presently shown), of which the former only is appropriate to the case of a 
solid sphere, the other becoming infinite at the centre. We have therefore to choose 

T> — Q 

±.\j/i — k%. 

With regard to the integer n, it must be at least as great as m, but may have any 
value from m up to oo. The form of the solution being now ascertained, the values of 
X may be found from the condition that at the boundary 

dH 
R # =0 or -ry+^R^=0 when r—r , 

and the arbitrary constants A, B . . . may be found from the initial distribution. 

4. We shall now determine the appropriate transformations of (1) and (2) for an 
ellipsoid of revolution, and shall confine ourselves in the first instance to the case of an 
ovoid ellipsoid, reserving that of a planetary ellipsoid for subsequent treatment. 

The axis of z being that of revolution, put x=p cos </>, y^p sin </> ; equations (1) and 
(2) become 

jVIJ-)OOOij-K.-X A.. !\i 



l 2i"Zi 
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dY dY 
V=0 or l^r+n 1 : -+tjV=0, 

in which I n are here the direction cosines of the normal to a meridian section. Let 
us now replace p and z by 

p=G sinh a sin ft, 2= c cosh a cos /3 ...».*. (11) 

where a and /3 are the thermometric parameters of the confocal system which includes 
the principal elliptic section of the bounding surface for which a:=ot 0? and whose axes 
are therefore 2c sinh a , 2c cosh a . 

Since p and z are conjugate functions of a, f$ 9 






da) + W / 



+ ~ 9 7)=c v (cosh "a — cos 73) ( ^ + 



c£/) 3 (^ 3 / 



\ dp 2 dz 2 



Moreover 



whence 



dV . dY . c£V 

4-=c cosh a sin /3 - +e sinh a cos /3 ■■■• , 
da dp dz 



dY . t ~<r£V 1 . ^dY 

-77: ~c sinh a cos p— c cosh a sm p-T) 

d/3 d/> ? ~ 



& 



rfV . dTV . dV 

e(cosh 3 a sin 3 /3+sinh 3 a cos 3 /3)—= cosh a sin j8— + sinh a cos j8— , 

ttO u>5t u>p 



and 



cosh 3 a sin 3 /3+sinh 3 a cos 3 /J=cosh a a — cos 3 /3. 



By the help of these formulae we may transform the general equation of conduction 
into 



dt 



d*V , d?Y , xl dV , , n dV , / 1 , 1 \d*V c 3 1 Q m ^V 

+ 7-r + coth a— +cot p- ? - + TT^+'^Tr "T^^ i"(cosJr' a— .cos" p)— . (12) 
da 2 df& dot d/3 \sm 2 /3 sinh 3 a/ ^</> 3 t x ^ } <*+ x ] 



In a similar manner, the equation to be satisfied at the boundary becomes 



dV 



V=0 or else— -\-^c^/cosh 2 a—cos 3 /3V=0, when a=a . . . . (13) 



da 



We must also find the space element, dE, 



cZE = pd&.dpdz = pd<f)dad/3 



dp 


dp 


da 


d0 


dz 


dz 


da' 


d/3 
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and finally, 

dE=c s sinh a sin /3(cosh 3 a—cos?ft)d(j)dad/3 . . . . . (14) 

5. We proceed to find the solution of (12) which is appropriate. We may satisfy 
it by putting V=2(cos m<£IT l + sin m<£U 2 ), where U^ and U 2 both satisfy 

tt+t^ +coth a-, — hcot p~—m z [ "-rr~+ T77. U = rfcoslr a—cos 2 p)— . 
da 2 ^/3 2 «« ^/3 \sinh 3 a ! sm 2 /3/ t x ' 7 dt 

And with regard to m, it must be observed that it cannot be other than a whole 
number, since the value of V must repeat itself in going round the surface of the 
ellipsoid in the ^-direction ; that is to say, 

lit \j , -L , Zj • • • "J^ • » • • • • • » * • I It/ J 

We may also put 

(16) 



• *#••*••• 



v=^J(/3).Q m k (a) 

where 3- and XI are functions of /3 alone and of a alone respectively, determined from 
the equations 

f|+cot^-^?-^=:XVcos 3 ^--^ (18 y ) 

dp 2 dp snr p v y/ 

— o+cotha— rrr- 12= —X% 2 cosh 2 aO+ to (19) 

n^ 3 rfa Sinh 3 a ' V '' 

wherein k is a constant, as yet undetermined. In the sequel it will appear that k 
has an infinite number of values for a given value of m and a given value of X, and one 
of the objects of the present investigation is to furnish the equation which determines 
it, and to approximate to its different values when the eccentricity of the ellipsoid is 
small. In this respect the present problem differs essentially from the corresponding 
one for a sphere in which k is independent of X ; it is then given by 

k=n(n-{-l), where n=m, m+1, m+2, . . . oo (17) 

If we put cos/3 =i>, the equation in & may be written 

( 1 .^)^_2 l ;^_ 7 ^-A=X» c VA-fe9- (18) 

v ' dv* dv l—v* x ' 

and, if we write £=\c sinh a, the equation in O is 

R 2 
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(\*c*+g*)a+kn 



or, more conveniently, 



2 






+2f-^+f o~to+xv(5+ 1 ^° 



wr^ m 



+o 



m/ 



. . . (19) 



It will be observed that these transformations are suggested by the consideration that, 
when c=0, the equations in 3- and 12 should reproduce the corresponding equations in 

P OT * and B^ for a sphere ; - is the semi-axis minor of the confocal ellipse, and we shall 



X 



V 



use - to denote in like manner the semi-axis major. 

6. I shall now show how equation (18) may be satisfied by a series of associated 
functions of the order m. If we compare the functions P,„ w+1 , P,/, P*/~\ in which the 
constant multipliers are so chosen as to make the coefficient of the highest power of v in. 
their rational factors unity, it is known that 



p M_p «+i i 



n* 



lib -J-^y "I 

L) «— -J. 

4*i* -1 w ' 



and hence 



o-r^ ., -T-. j.o . 27l 3 + 27^ — 2m 3 -— 1^ , (m 2 — m 2 )(n — l 2 — m 2 ^ ~ 



(2ra-l)(2» + 3) 



(4w 2 -l)(4.w-l 2 -l)" 



» » 



. . (20) 



It appears from this that the last term will vanish both when n=m and when 
n=l-\-m, the theorem reducing then to 



Sp m = T> ffl+ U 



1 



2 m -f 3 



«?, 



f)-*- W7. J 



and to 



3 



7 .2p »>+l— p ^+3 i 

2m + 5 



p , 



w+l 



It is clear, generally, that z^ 2i P^/ can be expanded in a series of associated functions 
of an even or odd degree, according as n is even or odd. Bearing in mind that 



■(1-^-24-^+* 



civ 1 



civ 1 



ir 



P n - 

x m • 



-(n(n+i.)-k)V, 



11 



we can obviously satisfy (18) by the expression 
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fl/r/ P m a P m+ % 1 n P »»+4, 

or by the expression 



in which 



IW** 



m+l 



6,P 



m+3 



]-*- w 



i * * « 



I. 



l 






i'a«s=r(fi— *)«i-ao 

1 



I /v P «*+2rl \ 

* * • ITI'^f- 1 - m i * • * /? 



^XIU/^x f W -p » , . ^, 



II. 

1 



> 



• • ( 21 ) 



^ 



Pr + !<*'!•+ 1 = ""(^f — &) «;• — ^r-i P 's+ A+ 1 = X (*'s — &) &* — &*-i 



• • isUsLl} 



J 



wherein X s c 3 =€, 



2w 9 + 2% — 2m 2 — 1 , , . ..* (n 2 — m^)(n — i 3 — m 3 ) . ^ 

zc r =— 777- — ., X/0 , ox e+^(^ + 1), jP,.=77i — 1%//< - — ^ 9 — ~, w=m+2r 



(2rc-l)(2^ + 3) 



(4^2„ 1 )(4 % _ 1 2_ 1 y 






(2w-l)(2rc + 3) 



(4^_ 1 )( 4 ^_ 1 2_ 1 y 



The arbitrary constants C and D are introduced that there may be no loss of 
generality when we put one of the series a , a x . . . equal to unity, or one of the series 
6 , b x . . . . The two sets of formulae I. and II. above are precisely identical, making 
allowance for the difference in the values to be assigned for n, and therefore the con- 
clusions drawn from I. will be in general true also .for II. ; and on this account we 
shall confine ourselves more especially to the former of these two systems. 

We must first examine somewhat more closely the sense in which the series (21, a) 
satisfies (18): If we stop at the term CajP m m+2r in forming the expression 



(i_^_ 2 /^ 



WS 



dv' 



dv 



1 — v* 



+\w-m, 



it is not precisely zero, but equal to C(— ea t F m+2r+ ^+ea r+1 V m ^ 2i ) ; hence it is only when 
a r and a r+l are either zero or indefinitely small that the differential equation is 
satisfied. 

We have, therefore, to show that, for certain values of k which are definite in 
position, a r tends to zero as r becomes infinitely great. In other words, the values of 
Jc are the roots of the equation 

rr> '" * ' * * * * * * * * * V / 
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Similarly, when n—m is odd, the values of k are the roots of 

&„ = 0, (24) 

and we have to show that these roots converge to fixed values. 

But before entering on this discussion it will be convenient to take up the con- 
sideration of equation (19), and show that it can be satisfied by a definite series of 
known functions corresponding to the same values of h ; and, preparatory to doing so, 
I shall digress briefly into the solutions of the equation 



dm 2 dB ( 3 n(n + l)\^ . 



fj/y>" ty* fi/} /> 



My reason for doing so is that, for the transformations which follow, a connected 
view of the properties of these solutions is necessary. 
7. If we write x=\r, the equation may be written 

dor x dx \ xr j 

We shall now show that the equation is satisfied by 

E= &(-—) = S„, or by U=T u =x n - — - — .... (25) 

\x dx/ x J \x dx) x N 7 

If we write 'R=af.u m we find 



— r +2(n+l)— -^ — h^=0. 

If we differentiate this equation, we may readily put the result into the form 

d % /l du n \ , ^ / , ^ x 1 d (1 d \ , 1 ub M 

4 ! * l ^«+r-3T=°; 



+ 2 (»+ 2 )li( 



^ 3 \a? dx J x dx\x dx x dx 



and, comparing this with the former equation, we find 



1 d . x 

u *+v=-T» u « V 26 ) 



x dx 



We thus obtain, generally, w w =( - v J u Q . 
And it is easily seen that u is given by 
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or 



Clxr +^0^—0. 



+J3 3 



7 0' 



X 



X 



The two values of II are therefore those given by (25). 

We may readily find these values as expansions in powers of x 9 by putting t for x* 



i 1. sin \Z~t t cos \/t 

and expanding — y T ~ and 



t 



b. 



( — l) n X n 



H 



1 



1 



o 

X" 






1 



X* 



1.3 . . . (2t&+ 1) 1 l.(29i + 3) 2 l 1.2.(27i + 3)(2rc + 5) 4 



• » • 



T — ( ~ 1) " L3 ■ ' • ( 2n ~ *) 



14 



l 






7T + 



x' k 



l.(27i-l) 2 ' 1.2.(27i-1)(2ti-3) 4 



+ . 



•*> 



(27) 



Since the differential equation is unaltered by replacing nby -— (n+1), it follows 
that S.^! and T_ / ^ i are also solutions of it ; and in fact it is clear, in comparing 
corresponding terms, that 

k5—i»/_i = y I J L n 

the constants, introduced by integration in S_^ l3 being so adjusted as to make the 
first n-\-l terms agree ; no constants are to be introduced in determining T_^_ 1 . 

—fd\ n sin \/t 
The form S # =2^( — ] • — jj- also indicates the expansion as a Bessel's function, 

quoted by Lord Rayleigh, 



so 



k M — ( — l) w A/ if *J«+ir) 5 
T.=(-l)- Vf^W«>- 



The finite expansions of S„ and T„ are 

" — ( - 1.2 *2W + 1.2.3.4 '2V 



4 



,0? 



71 



1 nXn'-l'Xn'-V) 1 



1 2^ 



1.2.3 



3™4 



23/V>< 



+ . . . j COS(«!+— J 



2/ 



. / . W7r\ 

sm (#+— - ) 



^ 



y. (28) 



hn 



1 »'(»'- 10 1 



■n 



,X 



'n 



JL.iJ 



2™3 



2*x 






cos (^+ 



^7T 



i %v- 10(^-20 



1 2^ 



L «-J. O 



1 \ . /' 7l7r\ 

+ . . . sm x+ 



2-V 



2"; 



in which n stands for n(n + l), and n—r-=(n-— r)(n+r+l), 
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Lord Rayleigh has given another form for 8 Jt in terms of differential coefficients 
with regard to x only, but what we are principally concerned with here are the 
expressions (27) which show that S n alone is finite at the centre, and with equation 
(26) which gives rise to the following formula! of redaction. 

8. Replacing in (26) u n by S,„ we find successively 



€Uk5 y~si 11/ £NJ 



r/ 3 S 






X" 



If we substitute these expressions in the equation which S M satisfies, we find 



WaM 
hence 



and finally, 



2n 4- 3 

KJu + %1 ^U-\- 1 "T" &U ^* 



# 



2n 4- 1 



a? 



S^+S^i+S^^O, 



<zs 



( 2ft + ^if =(»+ ^S^-nS^ ; 



fix 



S 



« 



^»+2 



+ 



2S, 



■n 









x? (2» + l)(2» + 3) (2»-l)(2»+3) ' (2»-l)(2» + l) 
1 rfS« 7i + l ^ . S« . mS 



w/ Cm 



■X 



> 



(2?j + l)(2re + 3) 



+ 



«-2 



(2»-l)(2» + 3) ' (2»-l)(2» + l) 



_y 



(29) 



We have to substitute these values in 



vvtv if J \AhAj 



m 3> 



& p 



which may also be written 



n(n + T)— m 3 Q 1 dS n 

rynl ' rA pi <Y* 

iAj <aj IAhaj 



% 3 —m 3 



(% + l) 3 — -m 2 ~ 2# + 2% — 2m 3 — 1 Q .„ „„ ~ 



(2^ + 1) (2^ + 3) 



(2ra-l)(2tt+3) * ' (2^-l)(> + l) 



If we now write 



Q - 



X.O.O . . . \jLiib"~~"X.J 



((n — l) 3 -~ m 3 )(?i — 3 3 — - m 3 ) , . . (m-— l 3 —m 3 ) 






» * » 



(30) 






*U + 



2 + (2w-l)(2^ + 3) ^^(i^~3)(4^1 3 ---l) ^ l j 
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If we now turn to equation (19), and compare the result just obtained with (20), 
we see that we can satisfy the equation (19) by 

Where a , a v a 2 . . . are connected by the same relations as before — namely, I. of 
(22), and the values of h are still, as before, the roots of ^ = 0. 

It will be observed that this result depends on the identity of the forms of the 
right-hand members of (31) and (20). But the transformation (30) depends essentially 
on the hypothesis that n—mm even. For the values of h which depend on n— m odd, 
and give rise to the second class of expressions for 3-, this method fails completely ; 
in other words, distributions of heat which are not symmetrical about the equator of 
the ellipsoid cannot be represented by Xl-functions of the type we have just found. 
I was led however to expect, from other expressions which will be given presently, that 
the true form in this case was to be discovered by putting 



3 C 3 ~/ 



u*— — -r u* • « . * *•#•* \0 Lit 



After substitution and reductions, we obtain 



e^+^+en'-kn'+W^-y^+n'+^a')^. 



Now the value of the expression 



by (29) 



If we put 



(w + 2) a -to 2 q , 2n* + 2n-2m*-l a {n-iy-m? 

: (2»+l)(2»+3) b " +3i " (2»-l)(2» + 3) *"'"(2»-l)(2» + l) "" 



~ J..O.0 . . , \4u7b mmm A.) „ 



where n— m is an odd number, w T e find 



d? 1 d 1— m 2 \^v _^ 2# + 2%— 2m 3 — 1 ~ (% 3 — m 3 )(% — l 3 — m 3 ) ~ 

,rff ~| ^ + + "^"/ y,l ~ yil+s+ "(2^-l)(2w+3) y * + (^»-l)(4^=l 8 -l) ^* 

It appears therefore that XI' can be expanded in a series of the form 

MDCCCLXXX. S 
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9. But we may also express 12 as a function of 77, the major semi-axis of the confocal 
ellipsoid. And, if we turn to the original equation in XI (19,) and put therein 
r)=\c cosh a, we find 

w _ x v)fp+2,f-|^ n =-va+ta (190 



VI 



We shall first replace XI by W, where X2=(?7 2 — X 3 c 2 ) 3 W, the equation becoming after 
substitution and reduction, 

(rf- XV)™-y+ 2(m+ % y^+m(w+ 1)W= -t? 2 W+£W. 
If we now replace W by w where W=^" w u?, we find after reduction 

and we can satisfy this equation by a series of S-functions. 
Before doing. so, we calculate the value of the expression 

d^S n 2m d$ n m(m + l)~ _____ ~ 2.(m + l) dS n n(n + l) + m(m + l)~ 
—— -— -|- - - q w — ^ — - — - |- - & n 

0,7)* 7] at) r\ A 7j ar\ rj" 

(n + m+l) (n + m + 2) ~ 2n*-±-2n— 2m 2 — 1 ~ (a, — m) (^ — m — 1) ~ . , , 

~~ j2n + l)(2n + 3) ^ +3 ~ (2^-l)(2^+3) b " + (2^~1)(2^4-1) ^~ 2 > bj W 

Putting 

1-Q.5 • . • (Jn—l) jy (qa\ 

1.2.6 . . . (n + m) A ' 

whether n— m be even or odd, we obtain 



<$ 2m d m(m + 1)\jj _tt 2% 3 + 2%— 2m 3 — 1 ^ (% 3 — m 3 )(% — l 3 — m 3 ) ^ 

^ ^ ^7 i?* - /*~ Mw+2 ~~ (2^-l)(2^ + 3y M ^ + (4^-l)(4.^.1 3 -l) H ^ 2 

The equation in w may, therefore^ be satisfied by either of the forms 






and m 

tt=~ L w 



> 



i . 



. (35) 
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It was the consideration of these formulae, including as they do both the cases of 
n— rn even and odd, which suggested the transformation (32). 

10. We shall now sum up the results of the last two articles and shall replace, in 
doing so, the symbols 2, Q, H by S. Neglecting unnecessary constants, we obtain 

(1) n—m even, distribution symmetrical on opposite sides of the equator, 






**(b ) — " ^0^ Omn -L. Q ^]V*» + 2 I 



1.3 



2m + 3 



(XoOj 



(2m + 5) (2m + 7)^ M+4> 

1.3.5 



(2m + 7)(2m + 9)(2m + ll) 



^3^^+6*T* 



m 



n(„)-(v*-W?L S , 2 ^ +1 ) a S 4- 2a(m+1)(m+2) a S 

, 2 3 (m + l)(m + 2)(m + 3) „, 
^(2m+7)(2m + 9)(2m + ll) s ° ffl+6 ^ 



1 



> . 



-J 



(36) 



(2) n—m odd, distribution equal and opposite on opposite sides of the equator, 



o_ ft p m+1 7 } p w+8 I 7, p «+B_7) p a»+7_L 



3 



^(b ) fc ~ 1 ^0"»+l "" 9™ , R ^l^w+s + A; 



3.5 



2m + 5 



(2m + 7) (2m + 9) 



60S 



2^+5 



3.5.7 



2m + 5 



(2m + 7) (2m +9) 



"^ 



(2m+9)(2m + ll)(2m + 13) & 3 S -+7 + • • • 



y. (37) 



. 2 8 (m + l)(m + 2)(m + 3) , q 

+ (2m + 9)(2m + ll)(2m + 13) ^+7 + • * • 



It will be observed that all these expressions for n, notwithstanding the factors in 
the denominators, are necessarily finite at the centre. 

The first equation of last article (19 7 ) suggests yet another form for 12 which we shall 
find of service in calculating the coefficients. 

Putting £= cosh a =~" 5 that equation may be written 



,2 



and a comparison of this equation with the corresponding one (18) in 9- shows that we 
may satisfy it by 
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or by 



Xl=a P„»(O-%EV K+3 (£)+ • • • +(-l) , 'a,P» w+3, '(0+ • • • 



O h "P M+l 






> . . . (38) 



according as n—m is even or odd. 

The function P w *(£) maybe written (£ s — 1) 2 ^}/(£) and we know that 



(ti + m) ! 



=r2 j 



^ w V ' 2»ml 1.3.5... (2rc-l) 

(m + 1) (m -f 2) . . . (m + r) 
(2m + 2r+l)(2m + 2r + 3) . . . (2m + 4r-l) 



9 n=m-{-2r 



and 



= 2 r 



(m + l)(m4-2) . . . (m + r) 



(2m + 2r + 3)(2m + 2r + 5) . . . (2m + 4r + l) 



, n=m + 2r+l 



This relation enables us to find the constant factor by which U(£) differs from 12(f). 
For when £ is infinitely near to unity, 



m 



O(0=n^ ¥ («^-"(l)"«iP-" +B (l)+ • • • ). »-w even 



m 



=^^'(6oP.- +1 (l)-6iP-- +8 (l)+ • • • ), n-m odd, 



while 12(f) then reduces to 



1.3 . . . (2m + 1) 
\c& g-»(-l) w+1 



, n—wi even 



1.3 .. . (2m + 3) 



, n— m odd, 



terms above f * +1 being neglected in both sets of formulae. 
Hence if 

12(f) =K.12(£), w— m even 

= K'12(£), n— m odd, 



> 



. . « . 



. . (39) 



( — l)»\'V s a K=1.3 . . . (2m+l) 



2(m + l) , 2 2 (m + l)(m + 2) 

1 ^(2m + 5)(2m+7) 2 



a ° 2m + 3 



2 3 (m+l)(m+2)(m+3) 



(_l)«+i\»+i c « +1 6 K / =1.3 . . . (2m+3) 



-\ 



(2m + 7) (2m +9) (2m +11) ^"^ ' " 

, 2(m + l) 7 , 2 2 fm + l)(m + 2) . 



2m + 5 J ' (2m +7) (2m +9) 



2 



2 3 (m + l)( m + 2)(m + 3) , 
(2m + 9) (2m +11) (2m + 13) 8+ ' " " 



(40) 
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11. We now return to the consideration of equations (22), confining ourselves in the 
first instance to the first system which, for the sake of convenience we shall rewrite, 

i¥*i=7(Ko-%) 
1 

1 



1 

6 



in which e=X s c 2 , 



2^ + 2^~2m 3 -l , , , ,v (7i 2 -m s )(7i-l 2 -m 3 ) 

^—""T^ — ^77o — r^r € + n ( n + 1 ) 5 Pr = 77~* — 7777==^ — To w=m+2r. 
(2w — 1) (2n + 3) '*■/>* (4t& 3 — l)(4tt — l a — 1) 

We are to endeavour to discover the nature of the convergence of the series 
a , a l9 <%,... and the nature and position of the roots of the equation a oo = 0. A 
similar system of equations has been discussed by Heine in his 'Handbuch der 
Kugelfunctionen/ second ed., p. 406, and the following investigations are mainly 
modelled on the principles which he has used. 

(1) Treating the constant a throughout as positive, we observe that the series 
a a x . . . a r are all positive when &= — co , and alternatively positive and negative 
when Jc= + °° • Moreover, as in Sturm's functions, no change of sign is lost or gained 
by the passage of any of the intermediate members of the series through zero ; and 
since the whole series gain r changes of sign as h passes from -co to + °° , it follows 
that all the r roots of <x r =0 are real. 

(2) We shall now show that all the r roots of a f =0 lie below k t \ but, previous to 

doing so, we must inquire more closely into the values of p r and -(k t+1 — K r ). 



\<Y\r — / YYh ) ('Yb ■ 1 " 11% J 

(a) p r z=z J — =- — — -; hence p r lies below yq (neglecting the case of m=0) 

ano. jLj ^ ==0 o/-'f" — m ~\ q* 

and, as the present discussion turns upon having this quantity greater than 1 + i~ 6 - 
we must assure ourselves that this is the case. Now, since Xc is not always necessarily 
small, and n need not be large, this expression is not always > 1+yg- ; but by choosing 
r and therefore n large enough we can ensure that this is the case whatever value Xc 
may have, and the present proof will commence with such values of r as certainly give 
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1 ... 1 

-(K r _ ml —K r _%)> l+iV When this is true it follows, d fortiori, that — (#c,. —1 — kv_ 3 ), 

1 11 

-(kv_i — kv~4<) . • . ~(/c r — K fmml ), -(Kr+i — Kr) . . . are all greater than yf. 

Let us also suppose, for the sake of clearness, that r is even ; then, if we substitute 
k=K r -i in the equations, a , a l9 a 2 , . . . a f _ 3 are alternately positive and negative and 
each is less than yg- of the one after it ; we have, in fact, 

jPia T =— "(ic^! — ic )a 



and since when k=K r „ l , p r a r =—a r __ 29 a r will also be negative. But, in the same 
manner, when we substitute k=K r or any greater quantity, the series a a Y . . . a f are 
alternately positive and negative, and remain so as k changes from k t to + co . We 
infer, therefore, that all the roots of a ? — lie below k t and that one root lies between 

Kr-i an d K r . 

(c) Moreover, only one root of a f =0 lies between K rmml and K r . For when ^=/c f _ 1 
the series a . . . a r _ x have r— 1 changes of sign and retain these ever afterwards as k 
increases ; therefore the subsequent changes of k can introduce only one more change 
of sign into the series d , a x . . . a r . 

(d) When k=p, one of the roots of a,.=0, the expressions -(Kr+i—p) 9 -(Kr+2—p) . . * 

6 € 

1 

are all positive and greater than 1 +-r 6 -; and, for this value of k 9 p r+2 a r+2 =-(K r+l — p)a r+} , 

6 
1 

^ f+ 3^ +3 =-(/c r+2 — p)a f+3 — a r+1 . . .; thus a f+1 , a r+3 . . . have all the same signs and each 

is less than - x - 6 - of the one after it, and these signs are opposite to that of a r „ v But 
as k increases from one root p to the next p of a f =0, a r ^ 1 must have undergone one 
change of sign ; hence a f+1 , a r+2 . . . must have each undergone one change of sign. 
In other words, each of the equations a f+1 =0, a r+9 =Q . . . has one root between each 
pair of a r =0. 

(e) It is clear, therefore, that each of the functions a 2 , <x 3 . . . vanishes once for values 
of k lying in the intervals between — 00 , k 09 k x . . . and but once. Let us there- 
fore conceive these lengths cut off on the axis of k, and construct the curves a,. +1 ==/[(&), 
a r =f^{k), a r -i = f$(k) . . . When a r =0, a y+1 = — 16a r-1 ultimately when r becomes 

6 

very great ; and when a r+1 =0, a r and a y _ 1 have like signs and a r = — —ra r ^ v a r there- 
fore becomes indefinitely small compared to a r _ l9 and reference to equations (22) shows 
that then a r-1 is indefinitely small compared to a f _ 2 . . . . The points, therefore, in 
which the curve a r —f%(k) cuts the axis converge to fixed points as r becomes infinitely 
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great ; and, in the neighbourhood of these points, the functions . . . a,._ 3 , a r _ %y a r ^ l 
converge with great rapidity. 

The reality of the roots of a^ = is thus proved, and the convergency, for them, of 
the series of coefficients. Similar considerations will apply to the roots of b^ = and 
the corresponding series of ^-coefficients. 

We proceed to approximate to the values of k in series ascending by powers of e. 

12. Although the existence and reality of the roots of ^ = are thus established, 
it is a hopeless task to attempt to find them, generally. But if we suppose the 
ellipsoid of small eccentricity, and confine ourselves to those values of X which are not 
very great, we may treat XV or e as being a small quantity ; and, if we conceive the 
roots expanded in ascending powers of e, a few terms of the series will be sufficient. 

When c=0 we know from the corresponding solution in the case of the sphere that 
the values of k are given by k=n(n-\-l), where r=0, 1, 2, 3 . . . oo ; we are therefore 
to expect that the new roots will consist of series of the form 

where k Y . . . are numerical coefficients. And here it may perhaps be proper to antici- 
pate a difficulty which may occur to the reader. The roots, as given above, are 
expressed in powers of X, which are themselves as yet unknown, being determined by 
the conditions to be satisfied at the surface of the ellipsoid. And the values of X 
depend again upon the particular root selected. Thus an apparent indeterminateness 
presents itself, which however is only apparent ; and it will be seen, a little further 
on, that the roots form a perfectly regular series, and that we can always choose the 
pairs of values of X and k which correspond to each other. For the present therefore 
we shall assume X to be known, and proceed to calculate the various roots of the 
equation in k. 

If we write f instead of K r —k, and calculate successively a x -r-a , a 3 -f-a . . . 

we find 

€ 1 p 1 a 1 ^a =<^ 

<?PlP* a Z + <*() = 0O0l( 1 ~ €*— 






l-#lJ& r+ A 



€ 4 ' 



0o0i 0102/ J 



And, in general, 

e'PiPs . . . p f a f -r-a o =0 o 1 3 . . . 0f^ 1 (l--eViS 1 +e 4 V ->> iS 3 ; ---e 6 ^ 1 S 3 4- . . .) . . (42) 
Wherein 

Pi Pz Pr 



r Si is the sum of the quantities ,,,,,....,, 

0001 0102 0r-i<P' 

r S s the sum of the products of every two non-adjacent terms of this series, 
r B s the sum of the products of every three non-adjacent terms, and so on. 
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The effect of having only non-adjacent terms is that no <f> occurs in a higher degree 
than the first, as would otherwise be the case. 
The equation which determines k is therefore 



Ui • • • ^(W.Sx+^.Ss- . . .)=o 



(43) 



If we neglect e altogether, the several factors reproduce simply the values of k given 
by k=n(n-{'l), r=0, 1 . . . oo. 

Before proceeding to approximate more closely to these roots, I shall state two 
properties of these sums which will present themselves in the reductions, the proof of 
which is so simple that it is scarcely necessary to dwell on it — 



f hi — r~l^l + 



P, 



<£,_!</>,/ 



Pr q I q Pr-iPr 



r^z 



4>>'-l4>r 



<f>r-<i<f> l \-i<f>r 



These properties will be of great service in proceeding to a second approximation, 
and in calculating the coefficients. 

Let us concentrate our attention on the (r+l)^ root, for which an approximate 
value is given by <£ r =0, viz., 



Aj=n(?i+l)+~-^ — _ _ e . . . 



(2n-l)(2n + 3) 



. . * » I :fc ~fc J 



The equation may be written in full 



l_ e 2(A +7 &_ + . . . + V' 



A^l 0102 



T~ i i ~T" • • • 



0,_10v <f>r<f>r+i 



4" € 



.4 



.PlJPs 



+ 



#L#4 



_<£o010303 <£ O 010304 



■+...+ 



# 



■ + 



JP»*+1 



i?i 



</>r-i<^ <MHi/\0O01 



-p . . . -p 



# 



r~\ 



^-2^-1 



+ 



it being observed that no two consecutive ps can occur in the products. 
Multiplying the equation by <}> r9 it becomes 



</y 



Pr . # 



1>1 



L +f ±1 +wfr+ 



_^,._! <f>, +1 \4>o4> 



+ €*(•■• .)+■•• 



0. 



We see from this that up to e 8 the value of <f> r will be 



. .. =0 



Pr i 1V+1 



\9 »-l 9 r+l, 



(45) 



Wherein <j>' r -.i and </>' r+1 represent what <£ /w and ^ +1 become when we substitute for 
/c its first approximate value (44). 
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To obtain a third approximation, we must substitute this value of h in the above 
equation (43). 

The third series of terms may be arranged 



A Pr ( Pi 



Pz 



£r-l\<M>l ' ^1^3 



"*Y" • • « " m y" 



Pr-2 



+ 



Pr+2 



<p r -$(pr-2 <pr+i<pr+% 



4- 
I • • • 






Pr-l i Pr+3 ■ 
• I T i i T • • • 



while the coefficient of (jf> f in the second series is 






,</>0</>l </>l$2 



<pr~2<pr-i 9^+ir*'+2 



+■ • • • , 



On putting, therefore, the value of <j> r given by (45) in these terms, we see that we 
shall have 

(46) 



L 21 P r I P r+l \ I M P r ~lP } ' i Pr+lPr+2 



9 r ~\ 9 y+i 



+ 



^',_<50 /2 ,_l fy'K+^'r 



•MM* 



+2/ 



where, in <£"r-i> <£"r+i we have to substitute the value of h given by the second approxi- 
mation. This gives h as far as e 5 . 

I proceed to find its value up to e 4 . 

Having n=m+2r, n=m+2r / , we shall put 

D f / }f =7?/(^+ 1) — -n(n-\- 1 ) 






2n' 2 + 2n f -2m?-l 2w? + 2n-2m? -\ 



'r'j 



(2w'-l) (2%' +3) (2t&-1)(2t& + 3) * 



and, for a first approximation, obtain 



7 / , ^\ , 2w 3 + 2w— 2m 3 — 1 
fc=n(ft+l)+ , 1wo „ , ON - e; 



for a second approximation, 
*=n(n+l)+ (2w _ 1X2w+3) - €• 



(2w-l)(2» + 3) 



J) ~*~D I * \T) 3 ' T) 3 



e 3 + . . . , (47) 



and, as a third, 



7 / , ,v , 2ra 3 + 2«-2m 2 -l 
&=«(%+ 1)H — — ztt^ — 7-zr- e- 



(2«-l)(2» + 3) 



+ e 



.4 



Pr 



P»'+i 



F* 






Pr+\ 









-\ 



Pt 



„ X-L'r— i,r Ur+i,r/ \*-' r—i,r -L' r+l,r/ -D >•—!,»• 



I IWl \ ff? O y-l,y ^- +1 Q y+1>| . 



D< 



♦•+V 



Pr—iPt 



Pr+iPr+% 



MDCCCLXXX. 



T) T) 3 D 3 T) 

T 



— ijiiiiiM 



(48) 
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This formula becomes simplified for the first of the series of roots corresponding to 
a given value of m y that, namely, for which r=Q; in the above expression p r =0, and 
the value of h is 

To reduce these expressions we observe that 

D, >v — (n — n) (V + n + 1) , 
D, +1 ,=2(2*i + 3), D r+Sl ,= 4(2w+5), 
D,._ lt ,= — 2(2ti— 1), D,_ 2i ,= — 4(2w— 3). 



- ,, . 2# + 2w — 2m 2 — 1 1 / . 

Jb urtnermore, since —z ~~r^ —=M 1 

(2n — l)(2n+o) 4 \ 



4m$ ~~ 1 



(2^-1) (2^ +3)/' 



8 — 2 (4m 2 - 1)(V -tt) (W + n + 1) 

;V ~(27^-l)(2^ + 3)(2^-l)(27^ / + 3) 5 

& — 4(2m + l)(2m-l) * _ 8(2w + l)(2m-l)(2tt + 5) 



(27i-l)(27i + 3)(2^ + 7) ? ' + *''' (2tt-l)(2n + 3)(2tt + 7X2rc + ll) 

_^ 4(2m+l)(2m-l ) * _ 8(2m + l)(2m-l)(2^-3) 



(2^-5)(2%-l)(2?i + 3)' '-*' (2w-9)(2ti-5)(2w-1)(2w + 3) 



Wherein, as before, w=m+2r. 

We may also write p in factors as follows :— 

(n— m )(n—m — l)(n + m)(n-\-m—l) 

Pr=Z (>-3)(2^--l) 2 (2^ + l) ' 

(w + 2 — m)(n + 1 — m)(^ 4- m + 2)(% + m+ 1) 

^ ,+1== "™ (2w + l)(2w + 3) 3 (2^ + 5) " 

Putting in these we find, for the first root of the series, 

1c —m[<m±X\±-±- 2(m + 1) € 2, 4(m + l)(2m + l) 3 ( . 



and, in general, 



7 / . .\ . 2^ 3 + 2n — 2m 3 — 1 7 „ . 7 o , /r ~\ 

£ =w ( w +l)+-— __— — € _^ € a+J s€ 3 + (50) 



(27&-l)(2*i + 3) 
where 
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Jr — 1 



^ 3 = (4m 3 — 1) 



(n-—m + 2)(n — m-\- l)(n-\-m-\-2)(n-t-m-\- 1) (n— m)(n—m — l)(n-\- m)(n + m — 1) 
(2ti + l)(2w +3)»(2^+6) (>-3)(2^-l) 3 (2?i4l0 

w-m + 2)(%~m+l)(% + m + 2)(^ + m + l) 
(2^-l)(2^ + l)(2^ + 3) 5 (27iT5)(2^ + 7) " 

(% — ■m)(n—m—' l)(n + m)(n-\-m — 1) 



> 



^ 



(2^~5)(2?i-3)(2?i~l) 5 (2^ + l)(2?i + 3)J 



>• 



13. The roots of the equation in. h being thus approximately found, I proceed to 



k 



show how to calculate the coefficients. When we consider that the series for # _ 
reduces to its first term when r=0, and to the term a r Vj l when Jc=n(n-\-l) and we 
suppose the ellipsoid to become a sphere, it is clear that, for the first root, a must be 
the leading coefficient, and that for the (r+l)* A root the coefficient a r must be the 
leading one of the series. Taking the general case first, we observe that, since none 

of the expressions for — , ~ . . . J:=I , — contains <f) r , it follows that a Q contains no terms 

Lvr\ Cvr\ Cvr\ VUr\ 

lower than €, a x none lower than €~ l . . . and a r ^ l none lower than e, it being under- 
stood that a r is finite and of the degree e°. It is also true, as we shall see, although it 
is not so evident at first sight, that a r+l is at least of the degree e, a r+2 of degree e 2 , 
and so on. We shall work out the first three coefficients on either side as far as e 3 ; 
though in the subsequent reductions the coefficients of e 3 in a r+1 and a r-1 are too 
complicated to be worked out fully in general. 
We have seen that 



PlPt • • • Pf^T = ^l • • • <^-l(l— r-A€ 3 + r-lS a .€*+ . . . ), 





a. 



PlP*P& • • • Pr+l*'' +l ~Z l = <f>Q • • • <M 1 — fSi.^ + rSjj.^ . . . ). 



a 



But we have also proved that 



Pt 



Q Q J_ — l 



Pr a Pr-lPr 






and that, up to e 5 , 



4>r-i<t>r (j> r -^r-l^r 



I a/ Pr i Pr+l \ ■ 4/ Vr-xPr , J^± L^+2 \ . 



hence, putting in this value in the expression for a r+l , 



f\ • • • Pr+l € n —rO • • • 9r-li € ^// + € A' I ^ JL/ ~>-l^l -T . . . f > 



T 2 
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and hence, dividing by the expression for — , we find 

a 



a, 



r+l ^ 






a r \(f> r+l (p r+1 (p r+ fy 



. . . . (51) 



In a similar manner 



a 



But 



and 









Pr-iPr 



r+ 



&*r-\<f>r 



Putting these values into the series on the right-hand side of (52) it becomes, up to e 5 , 



6 flr+iPr+2 

9 »-+i9**+2 



and finally 



Q>r+% 



: 2. 



a, 



and therefore 



'r 



a 



'*+&__*. 



9 »"+ir »'+2 



(53) 



Lvip 



= e 



9 »-+i-9 *'+2'r *'+s 



-^ 



The series for - JLJL , — , — 8 can be easily found ; for 

iXiif w,j* U>i}> 






and therefore, 









=^€ T^~ + e 2 T7 



<£ r „! 8 ^Vi^'r-s 



,a --— /y> /ri £■* 



VV/t 



9 ,.^9 ,._ 3 



Expanding these expressions in powers of e, we obtain 
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\A/f _1_ "1 • \A/f — — 



D 



p 3 & 



+ 



P< 



W'T— 1 • Cv^'— — 



&r+% • w?'"" 



T) T) 3 



D 3 



^+i i i* 



?*+2 j ^"»-4-l,r > i 



p r € O r _^ r Jp r € 

MM 



-L^*— 1»»' -L'r+i.r -*-A"+2»>' -L'r+i,'/ 
Pr 



T) 3 \ T) D 






Cv^ , „ - .p • VVJ' ' ' 



t^r + 3 T~ Ctj- — — 



w^_o T" (Jvp — — 



IrTjPv— 1 ^ 



2 



-L'f-l.f J-V— g,r 



Dr+is Ur+2,r Dr+$ t r 
Pt'Pr-iPr-% 



!^+ *^lU+ . . . 

-A-V+V M'+2,r/ / 

£*=+£^U+ . . . 

1J r-l,r u r-2,r/ / 
™-j- • • • 






JJ, r~l ) r J 



.€ 3 + . . . 



• (54) 



The determination of the ratios of the remaining coefficients is thus reduced to 
algebraical manipulation : the last terms, however, of a rmml and a r+1 are very complicated, 
and I have not succeeded in giving them a simple form, but they may, of course, be 
found in any case where assigned numbers are given for m and n. 

For the others we obtain, on substitution, 



w^_J_1 • \A/f——' 



e (2m + l)(2m-l) 

2(> + 3) — (2«,-1)(2m+3) 8 (2m+7) 



e 3 +B 3 + _ . _ 



,2 



w^_|_Q • \Atp— — 



1. 



4m 3 — 1 



8(2w + 3)(2%+5) 2 (2?i-l)(2?i+3) 3 (2?i + 5)(2™ + ll) 



.€ 3 + . . . 



\A/f _1_ Q • wj' "~" — 



48(2»+3)(2»+5)(2»+7) 



+ • .. 



' / ^^=— ."] • W^» ■— ■■»■■. 



(% 3 — m 2 ) (^ — l 3 — m 3 ) 
(2^ + l)(2^-l) 8 (2^-3) 



4m 3 — 1 



2 (2^-5)(2^-l) 3 (2^-f3) 



^+C 8 c 8 + . . . 



t^J'^ ,0 • WJ' ' 



(% 3 — m 3 ) (% — l 3 — m 3 ) (% — 2 3 — w?) (n — 3 3 — m 3 ) 



8(2^ - 7) (2w - 5) 3 (2^ - 3) 3 (2^ - lf(2n + 1) 



-2. 



4. 



4m 3 — 1 
(2^"-9)(2^-l) 3 (2^ + 3) 



€ 3 + • • • 



vv* O • Ivj' """^ 



(w 2 — m 3 )(«,— l 3 — mF)(n— 2 3 — rnF)(ii— 3 2 — m 3 )(?i — 4 3 — m 3 )(?i — 5 2 — m 3 ) 
48(2% ^llj(2»^9J*(2» - 7) 3 (2?i - 5)*(2» - 3) 3 (2» - l) 8 (2w + 1) 



• t | • • • 



1 4. We proceed now to the numerical calculation of several of the roots and co- 
efficients. 

I. Let m = 0. 



(1) Let r=0, ?i=0 ; a is here the leading coefficient and we may put it =1. 
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1 2 4 182 

k= S e ~ 135 + 3l5/7 + 3OT £ ^ ' * * 

_e e 2 91 3 

ai— 6 189 + 2.3*.5 3 .7 3e + ' " ' 

€^ 1 

/v -— - ■ ^3 .. 1 

4. *jU a > O • O • 4. X 

a3= 43l5?7 €3+ ' • ' 

(2) Let r=l, n=2, a x the leading term put equal to unity. 

L ...L 1 ■ 94 g_ 21388 31 

k "i~2l e 3 3 .7 se 3 2 .7 5 ,ll e ' ' ' 
e 1 

a2== 14 + 3V7Ul + * • " 

__f!_ ■ 1 8 , 

as ""504" l "2.3.7 8 .9.15 e "*" ' ' ' 

2 4 

a ° 135 ^3 S .5.7 ^ ' * ' 

(3) Let r=2, w=4, leading term (% put equal to unity, 

1c— 90 J- - _|_ 77674 a 2805228 g 
~ + 77 e + 5.7 3 .ll 3 .13 e 7Ml 5 .13.15 ,e + ' • • 

1 1 

e 2 1 

jl j. jb'x ^ ( /til iiu.iy 

8 16 2 . 

ai= ~" i715 e — 3.5.7 s .ll € + ' " ' 



8 3 _ _ 32 

3 3 .5 3 .7 s * e 3 2 .5 3 7 5 .11' 



^O Q3K3>73 ,e R2K8>75 1i e I ' ' ' 



li. J_jetj JTi — x« 



(1) For r=0 ? w=l, leading term a put equal to unity. 

h — oil. i_^ i ._ -3 i 

n, — '5 5 3 7 5 5 7 9 " * * 



£ 1 „ . 31 

10 375^" ' 5 3 .7 2 .10.11' 



/y i n i _ , ______ j^iv I _, ^ ,^£$ I L 

^l 1A Q>7K C T K3f72lA11 C » * " ' 



- £. _ 3 a , 

a * — 70 2.5».7.13 ? "T" • • ' 
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(2) r=l, n=3, leading term a T put equal to unit} 7 . 

, _ _7_ 1064 3 _ 808976 a 

+ 15 € + 3*.5 8 .7.11^ 3^.5 5 .7.11.13 € + " " " 

a2= 18 "~5.3 5 .13 e + ' " ' 

_ _^_ _ 3 

a3 ~ 792 9 3 .10.11.15 6 "*"••• 

£_ _8_ 9 

a o — 7 5 s6 3 5 6 ,7 e ' ' ' 



(3) Letr=2, n=5 






III. Let m=2. 



(1) If 7*=0, 7i=2 9 leading term a put equal to unity. 

1 7 3.7 3 3.7M1 

ttl_ 14 7 8 .11 + * * • 

e 3 5 „ . 

3_ " 504 7M5.18 ^ ' * * 



(2) When r=l, »=4 



(3) When r=2, w=6 



^77 78.113.13 ^ * " * 



75 1534 

+ 165 6 "'" ll 3 .13.15.17 e3 ~ 1 " • 



15. The functions with which we have just been dealing belong to the first of the 
two classes of £. There is, however, a second class in which 3- is of the form 

The investigation of these functions will proceed on the same lines as those we have 
already treated, and the general formulae for the roots and coefficients will still hold 
true, with the modification that we must here put w=m+2s+l» where 5=0, 1, 2 . . . 

For a given value of m, the first root of each series is 

, ( , -w \_cy\j_ * § 6(m + l) 2 , 12(m + l)(2m--l) 3 , . 

*=(^+lK m + 2 ) + 2 m + 5 €— (2 m + 5)8(2m + 7) € + (2m + 5) 5 (2m + 7)(2m+9) € + '*' ^ 55 ' 



144 



PEOFESSOE 0. mvm ON THE CONDUCTION 



As the mode of formation of the coefficients has been already sufficiently illustrated, 
I shall confine myself to writing down a few of the smaller roots of this class. 



Whenm=0, s=0, h= 2+ - e-g-^-— cS+ . . 



■, 7 ,„ , 23 23114 . 

m=0, s=l, k=12+ — e+gT^Tii e + • • • 

« t ™ . 59 . 696718 „ , 
m=0, ,= 2, *= 60+^6+ €»+.. • 



3 



*»=l.«=0,*= 6+ ^ C-j^g 



4 8 



3.7 6 .11 



i i 7 o^ . S1 i 21192 a i 

m=l, 6—1, £=20+^ e + 373-nri 3 e + • • • 



_i _o l,,,, 27 , 73892 
m-1, s-2, A-42+ 55 €+- 3#5 ? 11Sil317 



e a + • • • 



«=2, S =0, £=12 + | €- § ^€»+g^ i 3^+ . . 

17 506 

wi = 2, s=l, &=30 + - € +3TiJ]^3 €2 + * • * 



The larger m and n become the more nearly will the first few terms of the series 
represent the value of the root. 

16. We shall now show how to find the types of heat-movement which take place 
when the surface of the ellipsoid is maintained at a constant temperature zero. The 
general equation of conduction is satisfied by any expression of the form 

V= (A cos m<j&+/3 sin m^)e^ l %^)^\a) ; 

and, in order that the temperature may be constantly zero at the surface, we must 
have 

JV'K, X, m)=0 (56) 

or, more fully, 



in which the leading term is (-— l) r a r B m+ ^ 9 where m-\-2r=n. 

For given values of m and n this equation has an infinite number of roots; the total 
number of values of X is therefore, apparently, triply infinite. Now, in the corre- 
sponding problem for the sphere, the equation is S»(Xr ) = 0, thef number being doubly 
infinite; and there is no theoretical difficulty in using this solution to approximate to 
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the roots of il w w =0, where the eccentricity of the ellipsoid is small. The precise 

manner in which each of the roots of the first problem resolves itself into several in 

the second is interesting. We observe that S^ +2f may be the leading term of 12, either 

as the first term of f2 » +2 r, the second of Xl 1 Wi+3r .. a , the third of 2 ^ +2r _ 4 , . . . We infer 

therefore that each root of S =0 or of S x =0 corresponds to one for the ellipsoid, each 

root of 8 2 =0 or S 3 =0 to two, each of S 4 = or S 5 =0 to three; and, in general, that 

each root of 

Sg»=0, or of S 2#+1 =Q, 

corresponds to (n+ 1) roots of the equation in X for the ellipsoid. 

If, therefore, the total number of values of \ in the case of the sphere be nN, the 

total number for the ellipsoid is — - — N. 

All the roots of S»=0 may be found without difficulty by the general processes 
given by Professor Stokes (Camb. Phil. Trans., vol. ix.) and by Lord Rayleigh 
(Proc. Math. Soc, vol. v., p. 119); and, therefore, those of Q, m H can be expanded in 
powers of c 3 . We shall confine ourselves to determining those of O °=0. 

Employing the expression for 12(77) an d putting a =l, we have 

^o*i 1^ 1 &"• irir 2 4<i" * * • "~~" ^ • • • • • • • • v*^' / 

If we neglect e this equation reduces to S (\a) = 0, whose roots are given by \a=^7r. 
Let the full value of Xa be £*7r+^e 2 +^+ . . . , e being the eccentricity of the 

ellipsoid =-. The elements of the subsequent calculation are briefly as follows :— 



S (Xa) = v -- cos itt-\- ( -~- — -rr-5 ) cos in.e 41 -^ . . . 

S 9 (Xa)=—Tr-^ cos in— •(- -ir-Ahe 2 COSW+ . . . 

I 10 105\ . 



8 #7T 4 A , 

35 * 525 ^ "* 

Substituting in the above equation, and putting equal to zero the coefficients of e^ 
and e 4 , we obtain finally 

I 777" 

Xa=w+o^ eS +77^(^ 3 +27)6 4 + (58) 
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17. When the ellipsoid cools by radiation, the equation to be fulfilled at the 
boundary is 

-+tjcV / cosh 3 a-— cos 2 /3V=0, when a=a . 



d 



or 



rfV , t) 



cl, 



;+ x Vl~ e 2 cos s /8V=Q (59) 



observing that — - — =-=6. 
° cosn & a 

In its present shape, the process of satisfying this condition is complicated. If, 

however, we neglect e 2 , the condition then becomes 

The appropriate form of solution in this case is 



V= (A cos m^+B sin 7ii(f)e mkHl 3- m u il m Jl 9 



X being given by 



' jl L +'A%/=0, when £=&. ........ (60) 



If, as before, we neglect e 3 , this equation becomes simply 

the same equation as found for the case of a sphere. It is, therefore, only when we do 
not neglect eV in the expression \/l— eV 3 , that we obtain results belonging specially 
to the ellipsoid. We must accordingly indicate how the problem in its more general 
form is to be dealt with. To do so we require various general properties of the 
^-functions, to the discussion of which we now proceed. It should be understood 
that these properties are of a purely mathematical character, and have nothing to do 
with the special series of values which the physical conditions of the problem may 
ascribe to X. We have seen that for given values of X and m, the different values of h 
are perfectly definite, and it may facilitate the apprehension of these values to bear in 
mind the approximations which have been given for them in powers of c. 
The following theorems are well known :— 

+1 t-i r* ./ 7 r, f +l /T-> „\o 7 ^ (n—m)l (n + m)l »• ' 

_ P.-P.-efc^O, Ij^fdv^^-^j - > ( L-!)V -^ *V . . . (61) 
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and if we denote 
following : — 



r+i 



-l 



&m&m'dv by (n 9 ri, m, A.), or shortly by (n 9 n), we obtain the 



\1l 9 U ) — CIqCI ojn/ + (i\ CI ij m ' -{-Ct%Q;cyJ m H -f- . . . 

(n, n) =(a )^+(a 1 )^ w+2 + . . . 



It may be easily shown that 



for the differential equations satisfied by S- m n y 3- m n ' are 



>■:... (62) 



(^ n') = 0, n not = ft' . . . 



. . . . (63) 



and therefore 






r n On 

a\ Lis r^f 



' + 1 



n 



2v 



9, 






m 



2 



-\ 



1-z/ 3 



A*=\ a cV* w — ife M #' 



^ 



m* 



^ 



<:fe 1 — z^ 9 



;y=XWy-i^' 



-> 



• • • • 



( A/^ — ■— A/^/ ) J fir kt QjV' 

• -1 



(l-v 2 )(^ M 



49* 



■# 






+1 



-1 



0. 



(64) 



From this result we infer that any function may be expanded in a series of ^-func- 
tions ; and, first of all, we may invert the series which expresses the ^--functions of 
the a-group. This is manifest algebraically ; for if we solve the equations for &° . . . & r 
in terms of T? m m . . . P w ** +2r , we obtain for any one of the P's the following : — 



where 



±m l -—Jo& ~Tjl& ~\~JZ& "4" • • ' ? 



fs{s> S ) 



■ + 1. n . > 

-1 



• » 



• (65) 



.V 



Now we may assume that any function of v (at least any series of powers of v com- 
mencing with v m ) may be expressed in one or other of the three forms SAP,/* +2r , 

%BP m m+2r+l , t(AP m m+2r +BP w w+3r+1 ), and therefore in ^-functions of the a-group, &-group, 

cPS- 
or a combination of these. In particular v 2 &, zA£, yj . . . will give rise to ^-functions of 



civ* 



clS 



the same type as S-, and v&, — ... to functions of the opposite type , In general, if 



dv 



F(v) = XCi3-\(ii)Ci=\ +l ¥(v)^dv ....... (66) 



The expansion of v^& may be obtained without difficulty ; for if we combine the first 
of equations (64) with 

tj 2 
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(W) 



Cv X. ™ 



n 



m 



dv* 



2v 



rZP 

tt'.i. « 



w 



m* 



dv 






n(n+l)V 



n 
m •> 



we obtain 



\ 2 c 8 [ + Wp^=(^-n(w+l)) f >P»dk 



fe-w(w+l))(-l) r a/> r , w=w+2r, 



in which p denotes any even integer. 
But 

p' — n p'T> m -/V p'T> ?w+2_l_ 

hence, remembering equation (63), 



= -{p,jp'}, say. 



(67) 



The expansion of z^/ therefore becomes 



Z^ fOr^ 



l 

e 



JO, Of M ^(l, lf"'^- 



(68) 



1 8. These investigations and developments place vis in a position to deal with the 
boundary condition due to radiation; we may always put 



ly/l—<W3-*=g *9-°+g 1 n 9- 1 +g 2 "9*+ • (69) 



and all the g's may be expanded in ascending powers of e % , beginning with e% except 
g n n 9 which commences with e°. The general equation of conduction is to be satisfied 
in this case by the series 

V=(A cos mcf>+B sin m^)6- A ^(C ^°f2 +C 1 ^ 1 O 1 +C^ 2 a 2 + ...),.' . {10) 

in which the same X occurs throughout. 
At the boundary, where £=Ab, 



'da, 



Co fe 0+ ^ 0o oj+0^o^i+C 2 ^2+ . . . =0, 
C ^°O +C 1 (^H^ 1 « 1 )+C^ 1 ^ 3 + . . . 



f d£l { 



*s 



=o, 



C < g i °n +Ctf i 1 a 1 +CJ-^+g i *n i )+ . . . =o, &c. 
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These are the equations which determine the different values of X and the corres- 
ponding ratios C : G l : C 2 : . . . . When we neglect powers of e 2 beyond e 4 , the equation 
determining X is 



00 00 



^ffAm, t n V^ Y rnofc ~* ; • • • • • (n) 

If the first power of e 2 alone is retained, this equation breaks up into 

^+9r r &r=0,^+gfn,= 0,&c . . (72) 

The solutions of these equations may be easily found from the corresponding results 
for the sphere, and the quantities g/ * . . have been already found. 

19. We shall not pursue this inquiry further, but shall now show how the arbitrary 
constants introduced into the solution may be determined. As already explained in 
Art. 2 of this paper, when the solution is represented by Y=XAe~~ Knt v, A is to be 
found from 

A^v*dE=^vY dE, 
and what we have to find is, therefore, 



v 2 sinh a cos /3 (cosh 2 a — cos 2 /3)d<f>d/3da (73) 

-]y o 



'Zirf + lrvo 

o 



If the surface is kept at a constant temperature zero, or if we adopt the simple law 
of radiation — +~V=0, we may take 

dg A 



v= cos m<f>S- m n £l 



n 

m 



as the type of solution : the above integral may then be resolved into the components 

3- . &dv, v 9j & . 3-dv, sinh aO . Slda, sinh a cosh 2 all . £lda ; and of these the two 
-l J-i Jo 'Jo 

former have been already found. 

For the accurate solution of the problem of radiation we must take for v an expres- 
sion of the form 

cosm^(c ( y9-°fl +c 1 * 1 O 1 +c a ^O s + . . . ), 

in which c^ . . . are known constants, one of which is to be put equal to unity. 
When we substitute this expression for v we perceive that we have still another 
integral to evaluate, namely, 

sinh all . Xl'aa. 
o 
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J 



To effect the integrations it will be most convenient to suppose XI expanded in terms 
of £ in the form 

But, as is well known, 



i 



w-'^^^^iiTvMTv-v-* • • (74) 



and may be readily evaluated ; the result we shall denote by j(n 9 ri). 
We thus find 



XV^t-LV^/ 







n u n*'d£= tt(—) r + s (a,a\+a' r a s )j(m+2r, m+ 2s) + %a r .a r j(m+2r, m+2r) 9 

r not =s 



(75) 



11* I2*c?£ may be found by putting n=n\ a r =a r . . . 
For the remaining integral, combine the equations 



(i-D 



<pp. 



m 



d? 



n 



dV » 



m 



m" 



&% i-r 



-± m — AC t, jl m H\7l"-f m 1 ) x. m 3 



whence 



XV 



i 

~e 





CC\L CvX n, 



?np n rd£=(k-n.n+i)^ aPj'^+|(i-C 2 )^/^-n~M q . . (76) 



JtSUt 12(^4) — ®j{y- m ""^y^-m \ • « « • 



1 



Taking these together, we arrive at the value of £ 2 O.D<i£. 

20. The preceding investigations have related to the laws which govern the move- 
ment of heat in an ovoid ellipsoid. For a planetary ellipsoid, we must take 

p=c cosh a sin /3, z=c sinh a cos /3 ....... (77) 



The equation which V satisfies is 



a*v . ^ 3 v 



dV 



,*+*+ - -s+ - ^+fe-s^= *<«* »+ **" «>i> 



i y 2 v c 2 



cZV 



to satisfy which we must put 



Y=cosm$e-^%/(l3)n a n (K) . . . 



f 9 



. . . (78) 
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where & and XI are determined by the equations 

!&>+ cot/3j^~~-7-j^=~W l cos 2 f33---k$ (79) 

dp* dp sm 2 p v ' 

d?£l , . . dQ, , m 3 ft x o o • i > ^ . i 

rr+ tanli a-- -4 tt~= — -aV smlr aft+tfft. 

&a* da coslr a 

The latter of these equations, on putting £=Xc cosh a, becomes converted into 

/>0 . xO o\^ 3 ^ , ^>^ . ^ 3 ^ , ^ 2 A- 2 C 3 ft />0 x OX , 7 / X 

(p+X 3 c*)^- 8 +2£— +y y + -y 2 -=-(^-XV-)fl+to. . . (80) 

These two equations (79) and (80) are of precisely the same form, as (18 y ) and (19) 
which we have been discussing, and differ from them only in having — -X 2 c 2 instead of 
+X 2 c 3 . The subsequent formulae will, therefore, be the same almost in every instance, 
and the investigation need not therefore be repeated. 



